Partial scaling transform of multiqubit states as a 
criterion of separability 



C. Lupo^ V. Man'ko^ G. Marmo^ and E. C. G. Sudarshan^ 

^ Dipartimento di Fisica deU'lIniversita di Napoli Fcdcrico II and Istituto Nazionale 
di Fisica Nucleare (INFN) Sezione di Napoli, Complesso Universitario di Monte 
Sant'Angelo, via Cinthia, Napoli, 1-80126, Italy 

^ P.N.Lebedev Physical Institute, Leninskii Prospect 53, Moscow 119991, Russia 
^ University of Texas, Center for Statistical Mechanics, 1 University Station C1609, 
Austin, TX, USA 

E-mail: lupoOna. inf n. it, mankoSsci . lebedev.ru, marmoSna. inf n. it, 
sudarshanOphysics . utexas . edu 

Abstract. 

The partial scaling transform of the density matrix for multiqubit states is 
introduced to detect entanglement of the quantum state. The transform contains 
partial transposition as a special case. The scaling transform corresponds to partial 
time scaling of subsystem (or partial Planck's constant scaling) which was used 
to formulate recently separability criterion for continous variables. A measure of 
entanglement which is a generalization of negativity measure is introduced being based 
on tomographic probability description of spin states. 
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1. Introduction 

Entanglement is a purely quantum characteristic of composite system and it creates 
quantum correlations among the subsystems of the system. Till now there is no complete 
understanding and protocols to detect entanglement for generic multipartite system both 
for qubit (qufiit) case and for continous variables like photon quadrature components in 
the multimode electromagnetic field. There exists partial positive transpose criterion of 
separability ISEll^' this criterion gives necessary condition for separability of composed 
bi-partite system. In the case of two qubits and qubit-qutrit system this criterion is 
also a sufficient condition for separability which allows to detect entanglement without 
ambiguity. For all the other systems of qudits the positivity of partial transpose gives the 
necessary condition for separability only. The partial positive transpose (PPT) was also 
shown jl] to provide the necessary and sufficient condition of separability for a pure two- 
mode Gaussian state of photons. From the point of view of physics the partial transpose 
of subsystem density matrix means partial time reversal for the subsystem constituents. 
Obviously the change of time sign for all the system is an admissible operation in the 
sense that the density matrix obtained as result of time reversal corresponds to a state 
that can be realised in nature. But in case of entangled state change of sign of time 
for a subsystem constituents only is not an innocent operation and it can provide as 
result the new operator from the initial density operator of the system which is not 
an admissible density operator. The obtained operator can have negative eigenvalues 
which make it unappropriate to serve as a density operator. The transformations of 
density operator which provide as result a hermitian positive trace-class operator (which 
can serve as density operator) form the set of positive maps IH] or dynamical maps 
[ZIIHI. The transpose transform is an example of positive maps. The positive map of 
initial state density operator can be reahsed by time evolution of the system: it is just a 
dynamical map. Recently it was pointed out [OlEIj that there exist positive maps which 
are related to other transforms, e.g. scaling transforms of time, Planck's constant and 
other universal constants. In jlQ] the scaling of time (or scaling of Planck's constant) is 
applied to a subsystem of a composed system to formulate a new criterion of separability 
for multimode photon states. This criterion extends the criterion based on partial 
transpose operation and contains continous parameters describing the positive map of 
density operators. For entangled states the scaling provides as result operators with 
negative eigenvalues. Due to this the partial scaling transform can serve as necessary 
and sufficient condition of separability of Gaussian states of multimode electromagnetic 
field which extends result of |4| for bi-partite systems only. One can understand the 
time scaling transform as change of momentum of the particle p —>■ Xp. Obviously such 
change of momentum induced the transform of angular momentum (spin) also. 

Till now the partial scaling transform of time (or Planck's constant) was not used 
to detect entanglement in multiqubit (multiqudit) system. The aim of this letter is 
to introduce for the multiqubit system the map of density operator based on partial 
scaling transform of time for the constituents of subsystems and to use the transform 
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to distinguish the separable and the entangled states. We also apply the description of 
qubit states by means of probability distribution (tomograms) to consider a measure 
of entanglement suggested in \12\ . 

The letter is organised as follows. 

In Section El partial transpose of density matrix is discussed. In Section El 
notion of bipartite and tripartite entanglement is considered. In Section |31 a scaling 
transformations for quddits is introduced. Some examples (Werner [T^, GHZ J3| and 
W states [ini) for two and three qubits and qutrits are studied in Sectional In Section 
inland Section [71 conclusions and further topics are presented. 

2. Partial transpose 

The operation of taking the transpose of a matrix defines a positive map: 

T : p -> p > ^ > (1) 

it is well known that this map is positive but not completely positive: by the Peres- 
Horodecki [Hj criterion this gives a necessary condition for a state of a bipartite system 
to be separable. Physically this map is equivalent to time reversal: 

t — > -t (2) 

Suppose we have a bipartite spin system: 5* = 5*1 x 5*2, this yields the decomposition 
of the Hilbert space of the whole system as a tensor product Ti = Tii ® 0.2- A density 
state for the system S is described by a density matrix i?, the decomposition of the 
Hilbert space as tensor product induces the matrix representation R = Raa,b/3, where 
the latin and greek indices refer to first and the second system respectively. One defines 
the operation of partial transpose that consists in performing the transposition only on 
the indices of, say, the second system: 

T2 = / ® T : R = Raa,bl3 = Raf3,ba (3) 

This operation can be seen as a partial time reversal that acts on the second subsystem 
only: if the states present only classical correlations the map is always well defined and 
provides as result an admissible density operator of composite system. Otherwise, in 
presence of quantum entanglement the map may transform a density matrix of composite 
system into a matrix that is not positive. 

3. Bipartite and tripartite entanglement 

Positive but not completely positive maps are known to provide entanglement witnesses 
The usual approach is to take a bipartite quantum system with Hilbert space of 
states Hi ® 7^2 and apply a map like the following one: 

V2 = I®P (4) 

that acts trivially on the first subsystem and by a positive but not completely positive 
map P on the second one. 
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Partial transposition is a special example of such a procedure, in this case we take: 

T2 = / ® T (5) 

In general for a given P the criterion is only necessary for separability of states. But for 
the whole set of positive but not completely positive maps, the criterion becomes also 
sufficient. This is not a simple issue because it is not yet known how to characterize 
positive but not completely positive maps for a system that is not a qubit. Nevertheless 
it is known that for a two qubit and a qubit-qutrit system ^ and for pure Gaussian 
states jl] the criterion is also sufficient. 

Even for qubit system the situation becomes very different when one deals with 
multipartite systems. Suppose for simplicity to have a tripartite spin system 5*123 = 
5*1 X 5*2 X 5*3, in this case we must specify what kind of entanglement we mean. 

Indeed, one can consider bipartite entanglement for each one of the three possible 
bipartitions of the system, namely Si x 5*23, 5*2 x 5*31 and S3 x Su, otherwise one could be 
interested in tripartite entanglement: this yields to the classification of states discussed 
in [IE!. 

We say that a state is tri-separable if density matrix admits a decomposition as a 
convex sum of local states [13]: 

j j 

otherwise the state is entangled. 

The following relation holds between bi-separability and tri-separability: if a state 
is tri-separable it is also bi-separable for each choice of the bipartition, but a state can 
be bi-separable for certain bipartion but not tri-separable. It was shown in ^\ that 
there are also states that are bi-separable for all possible bipartition of the system but 
are tri-entangled. 

It is clear that a map like ^ cannot be used in order to obtain tri-entanglement 
witness because it is sensitive only to bipartite entanglement. In order to investigate 
tripartite separability one has to consider maps like the following: 

V = I®P®P' (7) 

where P and P' are positive but not completely positive maps acting on the second and 
the third subsystem respectively. 

It is clear that it cannot be sufficient to consider only transposition in order to 
obtain tripartite entanglement witness, this is simply because / T2 ® T3 = / ® T23 and 
one can obtain only a bipartite entanglement witness. Thus one can consider at least 
two different positive but not completely positive maps. 

Let us suppose that we have three families of positive but not completely positive 
maps acting on the subsystems Pj(Ai) for i = 1, 2, 3 depending on an array of parameters 
that we indicate with the short hand notation Aj, and consider the map acting on the 
whole system 

Vx,X,X, = Pl(Ai) ® P2(A2) ® P3(A3) (8) 
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Let us also suppose that for a certain value of the parameters, say Aj = 1, -Pi(O) = /j 
thus we can study tripartite entanglement by varing all the parameters or bipartite 
entanglement by fixing two of them to zero. 



4. Scaling transform 



In a previous paper JUI a family of maps was presented that can be useful in order to 
detect entanglement in multimode pure Gaussian states [H QU] of electromagnetic field. 
It is defined in a continous variables system as a rescaling of momentum: 



(9) 



I X — > X 
\ p — > \p 

G [—1, 1] it defines a semigroup of maps that are non canonical almost 



For A 

everywhere. The map Q is physically equivalent to a rescaling of time t — > At or 
Planck's constant h — > Xh. In the case of A = —1 it reduces to time reversal while 
for A = 1 is the identical map. The fact that © is in general non canonical play a 
foundamental role in order to detect whether a given state is entangled ^0]- It the 
present letter we propose a realization of this map in the case of spin systems, that is, 
those quantum systems with discrete variables. In the case of a qubit a generic (mixed) 
state is written as 



P=2 



l + z 
X + iy 



X — ly 
l-z 



One can define: 



P 



Txp = - 



l + z 
X + i\y 



iXy 
- z 



(10) 



Tx defines a scaling transformation that for A = 
A = — 1 reduces to the transposition or time reversal. 
The matrix representation of the map is 



1 is the identical map and for 



1 

i(l + A) 1(1 -A) 



1(1 -A) i(l + A) 



















1 



(12) 



that is, the matrix Ta transforms the vector 
Pii 

Pl2 
P21 
P22 

which corresponds to qubit density matrix 



(13) 



Pii 

P21 



Pl2 
P22 



(14) 
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into a new vector 



Pa 



Pii 

l+A^ _|_ 1-A^ 
^~Pl2 H 2~P21 

^~Pl2 H 2~P21 

P22 



(15) 



This new vector corresponds to new density matrix 



Pa 



l+A^ I 1-A 
^~Pl2 H 5~P21 



1-A , 1+A 
~^Pl2 H 2~P21 



P22 



(16) 



In order to determine if this map is completely positive one has to consider the 
auxiliary matrix: 



A) 












1(1 



1(1 



2V^ A) 





i(l + A) 



1 



(17) 



1 



1(1 + A) 

The map is completely positive if and only if the matrix Bx is positive [3 and it is 
so only for A = 1. In analogy to the fact that (jH)) is non canonical, Tx is neither unitary 
nor completely positive for almost all values of A. 

Notice that Tx can be written as a convex sum of the identity map and transposition: 



1 + A 



1-A. 



-T 



1^ 



2 2 

Being a sum of a completely positive map and transposition, it is clear that Tx cannot 
be useful in order to detect bound entanglement, i.e. Tx can detect entanglement only 
if T can. As a consequence bound entangled states cannot be detect by the map Tx. 
Equation ()18|) can be taken as a definition of the scaling transformation for a qudit 
system. 



5. Some examples 

In this section we are going to test the power of the scaling tranform to detect 
entanglement in bipartite and tripartite systems. 

5.1. Werner states for two qubits 

In this section we are going to test the capability of the maps Tx as sources of 
entanglement witnesses in the case of a two qubit system. In order to do that we 
apply the one parameter family of maps: 

rx = I®Tx (19) 

to the well known Werner states [12] . Werner states for two qubits are defined as follows: 



Wp=p|^)(^| + (l-p)- 



(20) 
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where is a maximally entangled (Bell) state: 
l^) = ^(|00) + |ll)) 



(21) 



and 14/4 is the maximally mixed state for a two qubits system. The density matrix 
expression for a Werner state is: 



- l+£ 

4 








2 
2 





4 














l-p 

4 





p 

L 2 








1+p 

4 



The matrix ()22|1 is known to be positive for all p G [— |, 1], separable for p G 
and entangled for p G (|, !]• 

Applying the 1 ® map one gets: 

i±£ 





(22) 



- -1 

3' 3J 



2 H^) 



l-p 

4 

2 



Wp-\ = TxWp 



with A G [-1,1]. 

The four eigenvalues of the matrices Wp-\ are: 

1 -pA 






P ( l+A 


2 V 2 







(¥) 


i-P 





4 





1+p 

4 



(23) 



ei(p; A) = e2(p; A) 
e3(p; A) 
e^ip; A) 



1 - 2p + pA 
4 

1 + 2p + pA 



(24) 
(25) 
(26) 



The map T\ detects entanglement in Werner states if some of this eigenvalues is strictly 
negative. 

A simple calculation shows that Tx detects negativity of matrix Wp-x only if 
p > (2 — A)^^ hence Tx gives a necessary and sufficient condition for separability since 
for A = —1 all the Werner states with p > 1/3 give at least one negative eigenvalue. 



5.2. Measure of entanglement 

In P2] the unitary spin tomogram of qubit state was introduced. It is defined as 
joint probability distribution w of spin projections mi, m2 depending on unitary matrix 
elements m G U{A) group. The tomogram of a state of two qubits with density matrix 
p reads 

Wp{mi,m2,u) = (mi,m2 | u^pu \ mi, 7722). (27) 

The tomogram completely determines the density matrix p ^ Ej- The measure of 
entanglement introduced in ^21 reads 

A^(p) = max^ {\wp^,j{mi,m2,u)\ - Wpj^j{mi,m2,u)y (28) 

mim2 
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Figure 1. on the left: phase diagram for matrices H23(l : on the vertical axis the 
parameter p definig the state H2()|l is plotted, on the horizontal axis the parameter A 
definig the map ((T^ is plotted. The "sep" indicates the domain of separability of 
()22|l . the coloured regions corresponds to negative eigenvalues in H2^^|l that witness 
entanglement. On the right: the anologous phase diagram for two qutrits state H56(l . 



The function M.{p) = for separable states and it is positive for entangled states. The 
maximum is taken with respect to all unitary group elements and with respect to all 
positive but not completely positive maps M of subsystem density matrix which indices 
the map of composed system density matrix p — > pu- 

If one applies all the maps one will have for some of them the matrix pu to be 
nonpositive for entangled state p. 

Below we applied the maps which depend on parameters A. In view of this, the 
maximum with respect to unitary matrix elements can be achieved if one takes the 
eigenvalues of the matrix pu- The maximum with respect to all positive but not 
completely positive maps M of the subsystem density matrix becomes the maximum 
with respect to A. 

Below we consider the eigenvalues of the matrix pu and apply the constructed 
measure in the reduced form. 

For the trace preservation property the sum of the eigenvalues equals one for all 
values of A: 

^e,(p;A) = l (29) 

We can introduce the quantity: 

m(p,A) = ^|e,(p,A)|-l (30) 

that is zero when all the eigenvalues are positive and is strictly positive when at least 
one eigenvalue becomes negative. If we take the maximum over all value of A G [—1, 1] 

M(p) = maa;Ae[-i,i]{m(p, A)} (31) 

this defines a measure of entanglement. 

For the previous example this maximum is reached at A = —1 for all values of p 
and yields: 
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Here it coincides with negativity measure jH]. 
5. 3. States of three qubits 

In this section we investigate the more interesting case of tripartite system, that is, a 
system of three qubits. 

For a composite system of n qu-(iits one can consider the states: 



Wp = p\w'){w'\ + {l-p)-^ 
where \w') is a maximally entangled Greenberger-Horne-Zeilinger like state: 



\w') 



\i)2 (S) ■■■\i)n 



(33) 



(34) 



The states ((221) are known [19] to be entangled for p > p^nt and separable for p < p^nt, 
with 

1 

For a system of three qubits (j34j) is the GHZ state: 
1 



(35) 



\GHZ) = ^i\000) + \lll)) 
this yields the following density matrix: 



(36) 
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p _l_ 1-p 



(37) 



Because of the particular symmetric form of the states (jHHj) they are bi-entangled if 
and only if tri-entangled. In order to appreciate separately bipartite entanglement and 
tripartite entanglement we modify states ()33|) introducing an addictional parameter: 



<e = P\^e){^e\ + il-p) 



for 9 e [0, |], where 



V2 



(|000) + cose|lll) + sin^|110)) 



(3^ 



(39) 



These states have no more the symmetric form of ()33|) and we expect a different 
behaviour when describing tripartite and bipartite entanglement with respect to different 
bipartition of the system. 

The pure state |\E'e)) (\E'6)| contributes to the entanglement of and the identity 
tends to remove purity as well entanglement while p is decreasing from 1 to 0. 
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The state vector in (j39p is maximally entangled for 6 = for all possible choices of 
partitions of the system, while for = 7r/2 it is separable for the decompositions 5*12 x 5*3 
while it is still maximally entangled for the decomposition 5*1 x 5*23 and 5*2 x S^i. 

Notice that the states Wp g in ()38p are symmetric for the interchange of the first and 
the second system: this implies that we can consider maps with two parameters only: 

I®Tx®T^ (40) 

for A = 1 we can study separability with respect to the bipartition 5*12 x 5*3, for yU = 1 
the bipartition S31 x 5*2 (that is the same as Si x S23) and by varing both A and fi we 
can study tripartite entanglement. 
We consider the matrices: 

<e;A,M=a®^A®T>;, (41) 

w'p Q is surely entangled if there is at least one value of A and fi for which w'p Q.y^ ^^ is 
strictly non-positive. 

More precisely if there is some negative eigenvalue on the line A = 1 this witness 
a bipartite entanglement with respect to the decomposition 5*12 x ^3, if a negative 
eigenvalue is found on the line = 1 this witness a bipartite entanglement with 
respect to the decomposition S'31 x 5*2, finally a negative eigenvalue in generic point 
witnesses tripartite entanglement, while the point A = /i = — 1 gives informations on 
the bipartition Si x S'23. 

In the symmetric case ^ = the eigenvalues, as functions of p, A and /x, are: 



ei = 


|+'7 + f(l + A)(l + ,) 


(42) 


€2 = 


^+'"'^ ^(1 + A)(1 + /.) 


(43) 


£3 = 


^ + |(l + A)(l-/i) 


(44) 


£4 = 


(1 + A)(1-^) 


(45) 


^5 = 


^ + f (l-A)(l + /i) 


(46) 


£6 = 


^-^(l-A)(l + /x) 


(47) 


ez = 


^ + f (l-A)(l-/i) 


(48) 


es = 


^:"-?(l-A)(l-.) 


(49) 



In order to measure entanglement we can consider the function 

8 

"^(p; A,/i) = 51 |ei(P>0; A,/i)| - 1 (50) 

i=l 

this function is zero when the matrix w'p.y^ ^ is positive and is strictly greater than zero 
when at least one eigenvalue is negative. 
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Figure 2. three qubit state (|38|) : some phase diagrams for 6* = for several values of 
p. On the horizontal and vertical axis the two parameters A and fi defining the map 
H4(J|I are plotted. Coloured regions indicate negative eingenvalues of (|41|l that witness 
entanglement of state For p < 1/5 the state is fully separable. 



A measure of entanglement is obtained by taking the maximum over A and /i: 

M{p) = maxx,f,{m{p] A, /i)} (51) 

it is easy to show that this maximum is reached at point A = — 1 and /i = — 1 and that 
m(p; —1, —1) = m{p] —1, 1) = m(p; 1, —1), this yields: 

That is, the amount of bipartite and tripartite entanglement is the same. 

The general case 6' 7^ can be worked by numerical analysis and some results are 
shown in figures ()3l4j) . 

For example for 9 = lA and p = 0.4 the state is surely tri-entangled and bi- 
entangled for the bipartitions S2 x 5*13 and 5*1 x 5*23, but it could be separable for the 
bipartition 5*3 x 5*12. 
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Figure 3. three qubit state H38|) : some phase diagrams for 9 = 0.7. On the axis the 
map parameters (|40|) A and fi are plotted. The states are entangled in any sense for 
p > 1/5. 



5.4- W states and qutrits 

In this section we study the states: 



p\W){W\ + ^Is (53) 

where \W) is the W state 

|iy) = i= (1001) + 1010) + I100)) (54) 

The amount of tripartite entanglement is equal to that of bipartite entanglement 
for all the possible bipartitions: as an example the tripartite entanglement measure for 
p = 0.6 is 

M{p = 0.6) = m(0.6, -1, -1) = m(0.6, -1, 1) = m(0.6, -1, 1) ^ 1.47 (55) 
As a simple example for a system of two qutrits we consider the GHZ-like state 

p|^)(^| + ^/9 (56) 
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Figure 4. three qubit state H38|) : some phase diagrams for = 1.4. On the axis the 
map parameters (|40|) A and /i are plotted. The state could be separable for p < 1/5. 
For p = 0.4 it is shown that the state is tri-entangled and bi-entangled for x S23 
and 5*2 x 5*31 but it could be bi-separable for S3 x 512. 



p=0.36 




Figure 5. three qubits state H53I) : some phase diagrams for several values of the state 
parameters p H53|) . On the axis the map parameters (|40|l A and /i are plotted. 
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where now 

IV') = ^(|00) + |11) + |22)) (57) 

As in fl2()|l the states in ()5fj|l refer to a bipartite system, hence we need only one 
parameter: see the right side of figure (|T)). For the state is fully separable for 
P < 1/4. 

6. Outlook 

We have seen that the proposed scaling transform cannot be useful to detect bound 
entanglement, nevertheless it is useful to study multipartite entanglement. In future 
publication we are going to consider a larger family of maps increasing the number of 
parameters. It is known that for qubits the set of positive maps is described by 12 
paramenters PH]. One could ask if it is possible to choice n < 12 paramenters that are 
sufficient in order to witness the entanglement of any non-separable states. 

Intuitively one can consider partial scaling of time transform with changing of 
"rotation angular velocity" which corresponds to changes 

for qubits. Analogous transform can be obtained for orbital momentum. 

For /ii,/U2,/U3 G [—1,1] this is a three parameter semigroup of transformations 
that corresponds to non isotropic contraction of the Bloch ball. These maps 
are positive for all the values of the parameters but are completely positive 
only for points inside a tetrahedron in the space of parameters with vertices 
(1, 1, 1); (1, -1, -1); (-1, 1, -1); (-1, -1, 1) [HI. 

7. Conclusions 

To summarize the results obtained, we can point out that we have suggested a new 
criterion of separability for multiqubit state. The criterion is based on partial scaling 
transform of time (or Planck's constant) which provides positive but not completely 
positive map of qubit density matrix. Using the partial scaling map of density matrix 
of composite multiqubit system, one can detect the entanglement. The partial scaling 
criterion is reduced to partial transpose criterion for particular scaling parameter equal 
-1. 

We used the tomographic measure of entanglement which for the value of scaling 
parameter equal —1 coincides with negativity measure. 

The Werner states for two qubit system and some families of states for tripartite 
system and bipartite qutrits system were investigated and the values of the state 
parameters for which the system state is entangled were determined using the suggested 
criterion. 
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In future publication, we study other examples of multiqubit states using the 
positive partial scaling transform criterion with increased number of parameters towards 
a complete family of separability criteria. 
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